First-order nature of the ferromagnetic phase transition in (La-Ca)Mn03 near optimal 

doping 



C. P. Adams, 1 ' 2 * J. W. Lynn, 1 ' 2 V. N. Smolyaninova, 2 ** A. Biswas, 2 R. L. Greene, 2 W. Rat cliff II, 3 S-W. Cheong, 3 

Y. M. Mukovskii, 4 and D. A. Shulyatev 4 
1 NIST Center for Neutron Research, National Institute of 
Standards and Technology, Gaithersburg, MD 20899-8562 
2 Center for Superconductivity Research, Department of Physics, 
University of Maryland, College Park, MD 20742 
3 Department of Physics and Astronomy, Rutgers University, 

Piscataway, NJ 08854 
4 Moscow Steel and Alloys Institute, Moscow 117936, Russia 
(February 2, 2008) 



Neutron scattering has been used to study the nature of the ferromagnetic transition in single 
crystals of Lao.7Cao.3Mn03 and Lao.sCao.2Mn03, and polycrystalline samples of Lao.67Cao.33MnC>3 
and La5/ 8 Ca3/gMn03 where the naturally occurring 16 can be replaced with the 18 isotope. 
Small angle neutron scattering on the x — 0.3 single crystal reveals a discontinuous change in the 
scattering at the Curie temperature for wave vectors below ~0.065 A -1 . Strong relaxation effects 
are observed for this domain scattering, for the magnetic order parameter, and for the quasielastic 
scattering, demonstrating that the transition is not continuous in nature, in good agreement with the 
temperature dependence of the central component of the magnetic fluctuation spectrum, the polaron 
correlations, and the spin stiffness reported previously. This behavior contrasts with the continuous 
behavior observed for the x = 0.2 crystal, which is well away from optimal doping. There is a large 
oxygen isotope effect observed for the Tc in the polycrystalline samples, and the Curie temperature 
is decreased by 7 K by substituting 50% ls O in the x = 0.33 sample. For the optimally doped 
x = 3/8 sample we observed T c ( 16 0) = 266.5 K and T c ( 18 0) = 261.5 K at 90% ls O substitution. 
Although Tc is decreased by 5 K for the x = 3/8 sample the temperature dependence of the spin- 
wave stiffness is found to be identical for the two samples. These results indicate that Tc is not 
solely determined by the magnetic subsystem, but instead the ferromagnetic phase is truncated by 
the formation of polarons which cause an abrupt transition to the paramagnetic, insulating state. 
Application of uniaxial stress in the x = 0.3 single crystal sharply enhances the polaron scattering 
at room temperature. Measurements of the phonon density-of-states show only modest differences 
above and below Tc and between the two different isotopic samples. 
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I. INTRODUCTION 

One of the simplifying features of conventional 
isotropic ferromagnets such as Fe, Co, Ni, and EuO is 
the isolation of the spin system from the underlying lat- 
tice. This made these types of materials ideal candidates 
to investigate the critical dynamics, scaling behavior, and 
nature of second-order phase transitions. This situation 
contrasts dramatically with the colossal magnetoresistive 
(CMR) oxides, exemplified by Lai-zCa^MnOs, where 
the electronic, lattice, and magnetic degrees of freedom 
are intimately intertwined. 1-8 This coupling leads to a 
number of interesting physical phenomena, including or- 
bital ordering, charge ordering, polarons, stripes, ferro- 
magnetic droplet formation at small 9,10 and large x, and 
the dramatic decrease in the electrical resistance in the 
ferromagnetic- metallic state that is the CMR effect. In- 
vestigations of these phenomena and the close relation- 
ship between CMR compounds, relaxor ferroelectrics, 
and high temperature superconductors has made this a 
field of intense study. An additional strong motivation 



comes from the possibility of technical applications of 
CMR materials in a new generation of magnetic sensors, 
read/write heads, and for a variety of new spin-sensitive 
electronics. 

In the CMR materials the transport and magnetic or- 
der arc connected through the ferromagnetic double ex- 
change mechanism 11 while the J aim- Teller effect locally 
distorts the lattice that surrounds a Mn 3+ ion, 1,2 ' 8,12-15 
coupling the electronic and magnetic degrees of freedom 
with the lattice. The resulting ferromagnetic transition 
was found to be dramatically different from conventional 
isotropic ferromagnets, and was interpreted as a dis- 
continuous transition between the ferromagnetic-metallic 
state and the paramagnetic-insulting state. 16 In partic- 
ular, the spin wave stiffness was found not to renormal- 
ize to zero at the Curie temperature as it should in a 
second-order transition. Instead, the spin wave inten- 
sities decreased rather than increasing, and the spectral 
weight of the magnetic fluctuations shifted to a quasielas- 
tic spin-diffusive peak in the fluctuation spectrum that 
rapidly developed as T — > Tc and peaked near Tc- 16-18 



Concomitant with this spin-diffusive component is the 
development of diffuse lattice (polaron) scattering, that 
is correlated in a way that is consistent with the forma- 
tion of zig-zag, CE-type structural stripes. The inten- 
sity for these polaron-polaron correlations also increases 
dramatically as T — > Tc in a manner very similar to 
the spin-diffusion component, and the temperature de- 
pendence for both of these closely follows the resistivity 
of the metal-insulator transition. 5 ' 18-25 It has recently 
been observed that these polarons form a glass above the 
ferromagnetic-metallic to paramagnetic-insulator transi- 
tion, which subsequently melts at a higher transition 
temperature to a fluid. 26 ' 27 

Considering the broad range of compositions, dopings, 
and crystal structures that comprise the CMR class of 
materials, it is not simply the observation of polarons and 
charge ordering that is significant, but the establishment 
of a clear relationship between the metal-insulator transi- 
tion and polaron formation that is crucial. We have con- 
tinued our neutron scattering investigations using single 
crystal samples of Lao.7Cao.3Mn03 and Lao.sCao^MnOs, 
and polycrystallinc samples of Lao.67Cao.33MnC>3 and 
La 5 / 8 Ca 3 / 8 Mn03 that have been treated with 16 and 
18 isotopes. The results demonstrate a ferromagnetic 
transition that is clearly first order in nature, as inter- 
preted in the original measurements on this system. 16 
The results also make it clear that the observed Tc is 
determined not by magnetic interactions but instead by 
polaron formation; the formation of polarons truncates 
the ferromagnetic-metallic state. We also show that the 
application of uniaxial stress enhances the polaron scat- 
tering. Finally, we have measured the phonon density- 
of-states both above and below Tc for both isotopes of 
oxygen, and observe small but measurable differences in 
the lattice dynamics. 

II. EXPERIMENTAL DETAILS 

The single crystal samples of Lao.7Cao.3Mn03 and 
Lao.8Cao.2Mn03 were grown by the floating zone tech- 
nique and have a mass of w 0.7 g. The crystal mosaics 
were single peaked with an intrinsic width less than a j 
degree. In this range of compositions the crystal struc- 
ture is orthorhombic but the distortion is small and the 
crystallographic domains appear to be equally populated. 
Therefore, on average, the sample itself is cubic, and we 
use cubic notation for simplicity. Nearest-neighbor man- 
ganese atoms are then along the [100]-type directions. 
The crystals were mounted in the (hkO) plane and the 
cubic lattice parameter was a = 3.867 A at room tem- 
perature. Elastic and inelastic neutron scattering mea- 
surements were collected on the BT-2 and BT-9 thermal 
triple axis spectrometers at the NIST Center for Neutron 
Research, using pyrolytic graphite monochromators, an- 
alyzers, and filters. Small angle neutron scattering data 
were obtained on the NG-7 30 m SANS spectrometer. 



The x = 0.3 single crystal has a Curie temperature 
of 251 K, which is the highest Tc for any Ca doped 
crystals. It has not proved possible yet to grow crys- 
tals with higher Ca content, and therefore bulk stud- 
ies at higher x are restricted to polycrystalline samples. 
For the isotope studies a 20 g polycrystalline sample of 
La5/ 8 Ca3/ 8 Mn03 was specifically prepared for the mea- 
surements, using the standard solid-state reaction tech- 
nique. This Ca concentration was chosen because it is at 
the maximum Tc and therefore any change in the Curie 
temperature cannot be caused by small changes in the 
overall oxygen content. 28,29 At each stage of the measure- 
ments high resolution powder diffraction measurements 
were carried out using the BT-1 powder diffractomcter 
to characterize the powder sample (s), and at each stage 
the samples were found to be single phase. To imple- 
ment the isotope substitution, the original sample was 
divided in half and each half was heat treated in par- 
allel with cither 16 O and 18 O, respectively, through an 
established procedure. 12 The fraction of isotopic substi- 
tution was determined by measuring the weight change, 
which established that approximately 90% substitution 
of 18 O was achieved. Bulk magnetization measurements 
were taken using a commercial SQUID magnetometer, 30 
while the ferromagnetic order parameter was measured 
by neutron diffraction. The nuclear scattering ampli- 
tudes for 16 O and 18 O are indistinguishable so that the 
total oxygen content can be determined by diffraction, 
but the fraction of substitution cannot be determined 
by this technique. After the 18 inelastic and diffraction 
measurements were completed, the sample was converted 
back to 16 O and the properties were re- measured to es- 
tablish that there was no change in the composition or 
changes in properties from the initial sample, within ex- 
perimental error. Finally, measurements were also car- 
ried out on the original x = 0.33 powder sample 16 ' 31 for 
comparison purposes. This sample was re-heat treated 
at higher temperatures compared to the original prepa- 
ration (and a small amount of Ca may have been lost 
in this process), and then the isotope substitution was 
performed. In these initial measurements only 50% sub- 
stitution was achieved due to the limited amount of 18 O 
available at the time, but the shift in Tc was still larger 
than at optimal doping. 

BT-2 and BT-9 thermal triple-axis spectrometers were 
employed to determine the magnetic order parameters 
and spin dynamics. For doping x between x ~ 0.15 and 
x w 0.5 Lai_a;Ca x Mn03 compounds are ferromagnetic at 
low temperatures. 28 ' 29 At smaller x there is anisotropy in 
the dispersion relations and additional excitations, 32 but 
in the regime of interest here the system behaves to a very 
good approximation as an ideal isotropic ferromagnct. In 
polycrystalline samples this makes it possible to measure 
the spin dynamics around the forward scattering direc- 
tion using the (000) peak reciprocal lattice point, 16 ' 33 
while of course in single crystals the dispersion relations 
and lincwidths can be measured around any reciprocal 
lattice point. 34 In either case the long wavelength excita- 



tions in a ferromagnet should have quadratic dispersion 

E sw = Dq 2 + A (1) 

with a temperature-dependent spin wave stiffness D 
and a possible spin wave gap A that would arise from 
anisotropy. The spin wave spectrum has been found to be 
gaplcss within the precision of neutron inelastic scatter- 
ing measurements, 5 which indicates that these materials 
are to a very good approximation rotationally isotropic. 

The phonon density-of-states was measured on the 
Filter Analyzer Neutron Spectrometer (FANS) and the 
time-of- flight Fermi Chopper Spectrometer (FCS). Un- 
certainties quoted for all the neutron measurements arc 
statistical and represent one standard deviation. 

III. MAGNETIC PROPERTIES FOR X=0.2 

In Fig. 1(a) we show a measurement of the ferromag- 
netic (110) Bragg peak in the vicinity of the Curie tem- 
perature for the x = 0.2 single crystal. Above the Curie 
temperature only the Bragg peak due to nuclear scat- 
tering is observed, while below Tc the magnetic contri- 
bution to the peak intensity is proportional to square of 
the order parameter (magnetization). This magnetic in- 
tensity is well described by a power law, except in the 
immediate vicinity of the Curie point where critical scat- 
tering is significant, and these data were excluded from 
the fit. The fit gives T c = 181.04 ± 0.14 K, and an 
(effective) exponent (3 — 0.324 ± 0.007 which is close to 
the prediction of the three-dimensional Hciscnbcrg model 
((3 = 0.36 5). 35 The errors represent statistical uncertain- 
ties only, and no account was taken of a possible spread of 
transition temperatures since the critical properties are 
not central to the present investigation. This Curie tem- 
perature is in good agreement with the values of 175 K, 36 
178 K, 37 and 185 K 32 reported previously for single crys- 
tals of the same nominal composition. 

The spin wave spectrum has been measured at small 
wave vectors, and the spin wave stiffness coefficient D 
is shown in Fig. 1(b) as a function of temperature. 
Each constant-g scan was least-squares fit to a quadratic 
dispersion law convoluted with the instrumental reso- 
lution. Correcting for the spectrometer resolution has 
the effect of decreasing the observed value of D by 
~ 5 — 10%, and the low temperature value was deter- 
mined to be 46 ± 2 meV-A 2 for this single crystal. This 
is in good agreement with the values previously reported 
for this composition. 32 ' 36 ' 37 As the Curie temperature is 
approached the spin wave intensities were found to in- 
crease due to the increase in the thermal population of 
spin waves, caused by the combined higher temperatures 
and lower D(T) values. In Fig. 1(b) we see that D ap- 
pears to renormalizc to zero at Tc, as expected for a 
conventional second-order ferromagnetic transition. We 
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FIG. 1. (a) Measurement of ferromagnetic order parameter 
at the (110) peak for a single crystal of Lao.sCao.2Mn03 near 
Tc- Data are well described by a power law with an exponent 
that agrees with a three-dimensional Heisenberg model, (b) 
Spin wave stiffness for the same sample. D renormalizes to 
zero as T — > Tc, as expected for a conventional ferromagnetic 
transition that is continuous in nature. 

did not find the development of a central quasielastic 
component to the magnetic fluctuation spectrum, in con- 
trast to the behavior found in the metallic regime. 16-18 ' 37 
For the x = 0.20 composition, Dai et al. 37 also did 
not find a central component, but they did suggest that 
D(Tc) was finite, which we do not see in the present 
measurements. 

The overall behavior changes as x increases to the 
optimal doping range for ferromagnetism. The excita- 
tions at long wavelengths are still isotropic spin waves 
with quadratic dispersion, but as mentioned in the in- 
troduction D fails to rcnormalize to zero at Tc- In- 
stead a quasielastic peak appears in the fluctuation spec- 
trum, indicative of the type of spin diffusion that occurs 
above Tc- This was interpreted as a coexistence between 
the low-T ferromagnetic-metallic phase, and the high-T 
paramagnetic-insulator phase, with the phase fractions 



changing with temperature. 16 ' 38 The spin diffusion com- 
ponent then grows at the expense of the spin wave inten- 
sities. Measurements of the magnetic correlation length 
as one approaches Tq from above show that this length 
does not diverge at the transition as would be expected 
for a continuous transition. 16,18,39 This agrees with mea- 
surements of the paramagnetic scattering, and the cor- 
relation length is only weakly temperature dependent at 
~15 A through the transition. These features cannot be 
explained in the context of a conventional second-order 
ferromagnetic phase transition, but can be understood if 
one assumes that long range ferromagnetic order is not 
the sole order parameter of the transition and that the 
coupling to the other order parameter forces the transi- 
tion to be first order. 40 



IV. MAGNETIC PROPERTIES NEAR OPTIMAL 
DOPING 

For a conventional ferromagnet above Tc the magnetic 
scattering is concentrated around the Bragg peaks with 
g-dependent intensity that follows the Ornstein-Zernikc 
form 



J a 



1 



(2) 



The correlation length £ (= 1/k) in real space is the pa- 
rameter that diverges at Tc for a second order transition. 
Such a divergence was not observed in previous polycrys- 
talline measurements, 16,39 or in more recent single crystal 
measurements made around Bragg peaks. 17 ' 18 However, 
SANS measurements on polycrystalline samples are re- 
stricted due to the large amount of small angle metal- 
lurgical scattering from the crystallites, and therefore we 
have pursued SANS measurements on the x = 0.3 
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FIG. 2. Integrated intensity of the small angle neutron 
scattering for a single crystal of Lao.7Cao.3Mn03. The scat- 
tering changes in an abrupt manner at Tc = 251 K. 



single crystal where high quality data could be ob- 
tained over a wide range of wave vectors and temper- 
ature, and where we have also been able to measure the 
time-dependence of the scattering. 

The results are quite different from conventional 
second-order behavior, as (now) expected. The Ornstein- 
Zernikc behavior in Eq. (2) is maintained for q ^ 
0.065 A -1 with a short and weakly temperature- 
dependent correlation length as found previously. 18 For 
smaller q, on the other hand, there is a discontinuous in- 
crease in the scattering at Tc that cannot be described in 
terms of the Ornstein-Zernike form or the common Porod 
form (log / oc log q). 41 This scattering is strongly peaked 
in the forward direction, and dominates the intensity on 
the (two-dimensional position-sensitive) detector. Fig- 
ure 2 shows the temperature variation of this scattering, 
recorded as the total counts on the detector. This corre- 
sponds to an integration of the scattered intensity from 
roughly q = 0.001 A" 1 to q = 0.015 A" 1 . The wave- 
length is equal to 10 A and the detector was 15 m from 
the sample for this measurement. The q dependence of 
the intensity is typical for scattering from domains and 
domain walls, 42 and the abrupt nature of its onset is con- 
sistent with the transition being discontinuous. 
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FIG. 3. Time dependence of the observed SANS scattering 
for three different wave vectors at 250 K. Data are scaled at 
t = and shifted vertically to the same asymptotic value, 
and are fit to a decaying exponential with a time constant 
r = 6.6 ± 0.3 minutes. 

One hallmark of first-order transitions is that, even if 
thermal hysteresis is too small to be observed, they are 
accompanied by relaxation effects. Figure 3 shows SANS 
measurements at three different wave vectors as a func- 
tion of time, after warming to 250 K. The time for the 
sample to equilibrate near Tc is clearly quite long. The 
data can be described by an exponential decay with re- 
laxation time t, and the solid curve is a fit of this form 
to the data. Different values of q gave similar values of 
r within statistical uncertainties, showing that all length 
scales measured in these SANS experiments exhibit sim- 
ilar relaxation times. Warming was performed in steps 
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FIG. 4. (a) Thermal irreversibility of the (100) magnetic 
Bragg peak on warming and cooling. The data were ob- 
tained by changing the temperature in 1 K steps every 2.5 
min. (b) Time dependence of the intensity of the (100) Bragg 
peak, both warming to the set point of 250 K, and cooling 
to the same set point. Data are fit to an exponential re- 
laxation. Note that the time constants are roughly identical 
and the asymptotic values are identical within experimental 
uncertainties, so there is no genuine hysteresis observable. 

of 5 K every 40 minutes for these measurements, and 
the relaxation time was found to be maximized at Tc- 

One of the difficulties in unambiguously identifying 
this as a first order transition has been that, although 
D does not renormalize to zero, the magnetic contribu- 
tion to the intensity of the Bragg peaks appears to be 
continuous. In particular, although there are some "ir- 
reversibilities" observed on warming and cooling, 16 no 
clear hysteresis has been identified. The availability of 
the x — 0.30 single crystal has allowed us to investigate 
the behavior of the magnetic Bragg peak in detail. In Fig. 
4(a) we plot the intensity of the (100) Bragg peak, and 
we see a clear difference on warming and cooling. This 
apparent hysteresis is due to relaxation effects, and in 
Fig. 4(b) we show the time dependence of the scattering, 
after the sample temperature arrives at the set point. 
The curves are fits to exponential relaxation functions, 
which describe the data quite well. The fitted r versus 




230 



270 



240 250 260 
Temperature (K) 

FIG. 5. Relaxation time r as a function of temperature for 
the (100) Bragg peak. The short time constant found away 
from the transition is typical of the time it takes for a sample 
to equilibrate, but in the vicinity of Tc the time constant 
increases dramatically, indicating that the transition is first 
order in nature. 

temperature is shown in Fig. 5, on warming and cool- 
ing. Although the statistical uncertainties are large, par- 
ticularly when the total change in scattering between two 
temperatures is small, we can clearly see that r is maxi- 
mized near Tc- 

A further result of the time-dependent measurements 
is that the asymptotic value of the order parameter at 
any particular temperature is identical on warming and 
cooling, as shown in Fig. 6. Thus there is no true ther- 
mal hysteresis in the asymptotic time region. We note 
that the value for Tc determined in these measurements, 
251 ± 3 K, is slightly lower than the preliminary value we 
reported previously 257 ± 1 K. 18 The (correct) value of 
251 K is determined here by locating the inflection point, 
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FIG. 6. Intensity of the (100) magnetic Bragg peak, in the 
long-time limit. The data reveal that there is no genuine hys- 
teresis observable. The bottom curve shows the derivative of 
the data, where the minimum identifies a Curie temperature 
of 251 K. 
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FIG. 7. Intensity of the quasielastic magnetic scattering for 
a temperature of 266 K at a wave vector of 0.11 A -1 . Note 
that it takes about an hour for the intensities to become equal 
on approaching from low versus high temperatures. In the 
long-time limit there is no genuine hysteresis observable. 

after taking the data on warming and cooling and prop- 
erly taking into account the long time-scale behavior of 
the system. This value is also in good agreement with 
the abrupt onset of SANS shown in Fig. 2. The SANS 
data also lack true thermal hysteresis. 

We also measured the relaxation of the central 
quasielastic component of the magnetic fluctuation spec- 
trum just below Tc for the x = 3/8 16 sample, and the 
data are shown in Fig. 7. The sample temperature was 
changed to 266 K from well below or well above this tem- 
perature, and the data collection at the elastic position 
was started when the sample block reached the set point. 
The measured relaxation times are somewhat different 
for the two cases, but it is clear that the asymptotic val- 
ues are identical within statistical uncertainties. Hence, 
the activation barrier and latent heat for this transition 
are rather small, which has hindered the identification 
of this as a first order transition in the specific heat un- 
til recently. 43 Previous measurements of the order pa- 
rameter (and measurements reported in the next section 
on polycrystalline samples) had always shown some ir- 
reversibility on warming and cooling as a result of the 
activation barrier, but the temperature variation in the 
order parameter always appeared to be continuous. 



A. Effect of uniaxial stress 

One of the characteristic features of CMR materials is 
their extraordinary sensitivity to external perturbations. 
The application of a magnetic field, for example, repre- 
sents a very modest energy in a system that orders at 
250 K, but a magnetic field can drive the metal-insulator 
transition and the associated CMR effect. Stress, which 
is always present in thin films, also has a dramatic effect 
on the electronic properties, and is one of the primary 



reasons for the rejuvenation of interest in this class of 
materials. 44 We have already established that the for- 
mation of polarons in the paramagnetic state is directly 
related to the resistivity, 18 and we may then anticipate 
that the polarons themselves will be directly affected by 
stress. 

To date there have been only a few studies of CMR as 
a function of applied force, using two-dimensional strain 
in thin films 45 and hydrostatic pressure. 36,46 In addition, 
the effects of internal pressure that result from substi- 
tution of smaller Pr ions for the La A-site cation, for 
example, have been performed. Other options for cation 
doping have shown that the behavior of the entire x — T 
La-Ca phase diagram can be achieved via changing the 
average radius at the A-site. 40,47 ' 48 Results can often be 
understood in terms of the double exchange model when 
the distortion of the ideal perovskite lattice is consid- 
ered. Changing the average radius of the A-site cation 
changes the tilt of the oxygen octahedra that surround 
the B-site cation, which effects the electronic band width 
and the double exchange. An increase in distortion away 
from the ideal perovskite structure tends to favor charge- 
ordered insulating behavior, where the antiferromagnetic 
interaction dominates over metallic ferromagnctism. 

Our single crystal of Lan.7Can.3Mn03 was mounted in 
a stress rig on the BT-9 triple axis spectrometer in the 
(hkO) zone. Stress was applied along the axis of the crys- 
tal, roughly 30 degrees from (001) in the direction of 
(220). All measurements were made at room temper- 
ature where the charge peak intensity has dropped to 
roughly 30% of its maximum value near Tc- Given the 
geometric limitations imposed by the stress rig we were 
essentially limited to measurements near the (220) peak 
using an incident energy of 52 meV. The charge order 
peaks are broad (corresponding to short range correla- 
tions) and occur at wave vectors (±^, ±j, 0). Near 
(220) these occur at (1.75 2.25 0) and (2.25 1.75 0); the 
structure factor is negligible for the "longitudinal" peaks 
at (1.75 1.75 0) and (2.25 2.25 0). In Fig. 8 we show 
a transverse scan of the charge order peak at (2.25 1.75 
0) at kbar and 1 kbar. The intensity of the peak has 
increased by nearly a factor of two, while the intensity 
of the fundamental Bragg peaks are unaffected at this 
level of stress. It is evident that the diffuse scattering, 
which is related to the "single-polaron" diffuse scatter- 
ing, has also increased. Uniaxial stress therefore has a 
rather dramatic effect on the polaron scattering. 

Because of the rather large change in intensity that was 
observed when 1 kbar of stress was applied, we reduced 
the force to zero in order to determine the behavior at 
lower values of the stress. The polaron intensity was re- 
duced, but not to the initial value; the observed scatter- 
ing was about half way between the initial data and the 
1 kbar data. Subsequent measurements up to 0.5 kbar 
showed no change, and with further increase the scatter- 
ing began to increase up to the 1 kbar data. Hence we 
believe that the scattering would increase approximately 
linearly with stress from the virgin state. Unfortunately, 
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FIG. 8. Change in the intensity of the charge-order peak 
when uniaxial stress is applied. The intensity increases by 
more than a factor of two in going from kbar to 1 kbar of 
uniaxial stress. The intensity of the single-polaron (diffuse) 
scattering also increases with applied stress. 

shortly after returning to 1 kbar of stress, the crystal 
shattered, ending the data collection. 



B. Oxygen Isotope Effects 

One very strong indication for some type of cou- 
pled transition in Lai-zCa^MnOs is the oxygen isotope 
effect. 12 ' 49-55 Conventionally, magnetism is a property 
that emerges from the electronic structure of the system, 
which is determined by the component elements and the 
crystal structure that they assume. Once that crystal 
structure is determined, the specific isotopic masses of 
the components do not usually enter into the electronic 
Hamiltonian. The classic example that runs contrary to 
this is BCS (Bardeen-Cooper-Schriffer) superconductors 
where the phonons lead to an attractive interaction be- 
tween the electrons and fundamentally change the nature 
of the low temperature electronic structure (a gapped su- 
perconductor rather than a Fermi liquid). Such phonon 
effects are limited to temperatures less than 40 K. How- 
ever, in Lao.8Cao.2Mn03 substitution of 18 for 16 leads 
to a 20 K decrease in T c , from 180 K to 160 K. 12 An iso- 
tope effect in antiferromagnetic Lao.sCao.sMnOa is also 
seen, but in the opposite direction and effecting the Neel 
temperature. 49 

To investigate the effect of oxygen isotope substitution 
on the magnetic behavior we chose to work at x = 3/8, 
where Tc is maximized and the CMR effect is large. 29 
This choice of x minimizes any possible variation in or 
distribution of Tc with overall oxygen content. Indeed, 
for the measurements in Ref. 12 the observed oxygen iso- 
tope effect has been partially attributed to different oxy- 
gen stoichiometry. 51 For the present samples the high res- 
olution powder diffraction measurements could be used 
to establish that the oxygen contents for both samples 
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FIG. 9. Measurement of the intensity of the mag- 
netic scattering on warming for poly crystalline samples of 
La5/ 8 Ca3/gMn03 that have been treated with ie O and ls O. 
The Curie temperature is reduced by 5 K for the 18 O sample. 

were identical within experimental uncertainties, and 
thus there can be no significant variation in Tc due to 
the total oxygen content. The lattice parameters, average 
atomic positions and occupancies were identical as well. 
Density measurements indicate that roughly 90% of the 
16 was replaced by 18 in the substituted sample. 

Figure 9 shows the intensity of the (100) powder peak 
as a function of temperature for the two samples, to de- 
termine the effect of oxygen isotope substitution on Tc- 
The data were obtained as a function of warming at iden- 
tical heating rates. The Curie temperature of the 18 
sample is 5 K lower than the 16 sample, as can been seen 
in the inset. The same 5 K difference was seen on cooling, 
along with a 5 K irreversibility in both the curves relative 
to the warming runs. Given the first-order nature of the 
transition the irreversibility is expected and is consistent 
with previous measurements. Taking the average of these 
warming and cooling runs gives Tc( 18 0) = 261.5 ± 0.5 K 
and Tc( 16 0) = 266.5 ± 0.5 K. As we pointed out earlier 
for the single crystal measurements, these Tb's are deter- 
mined from the inflection points of the order parameter 
curves. The value of Tc for the 16 O sample is indeed 
the maximum value seen on the Ca-concentration phase 
diagram and is a further indication of the high quality of 
these samples. 

This 5 K shift is much smaller than reported in the 
x = 0.2 sample, so for comparison we have also made 
measurements on the original x — 0.3 sample. With only 
50% 18 substitution we obtained AT C = -7 K. It is 
noteworthy that the transition for the partially substi- 
tuted sample was just as sharp as for the fully substi- 
tuted sample, indicating that there is no significant effect 
on the properties caused by the oxygen isotopic random- 
ness. The larger isotopic shift with decreasing x that we 
find is in good agreement with the trends already estab- 
lished for this system. 

Having confirmed that oxygen isotope substitution 
changes Tc, the next question is how the magnetic in- 
teraction strength changes, which we characterize by 
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FIG. 10. (a) Example of the spin wave dispersion relation 
for two typical temperatures. Data are well described by gap- 
less quadratic dispersion with the indicated values of D. (b) 
D versus temperature for the le O and ls O samples. The 
observed D values are identical at each temperature inde- 
pendent of the oxygen isotope, despite the difference in Tc 
(as indicated by the arrows). The extrapolated temperature 
where D(T) would renormalize to zero is 271 K, but the fer- 
romagnetic-metallic state abruptly disappears at the temper- 
ature where the polarons form, and this truncation occurs at 
a lower temperature for the 18 sample. The D(Tc) value for 
the original data in Ref. 16 with x = 0.33 is also indicated by 
the open diamond symbol. The inset shows the stiffness over 
a wider temperature range. 

measuring the spin dynamics. In Fig. 10(a) we show 
the \ow-q spin wave dispersion in the 16 sample at 
250 K and 265 K. The dispersion is well described by 
Eq. (1) with gapless excitations, typical of these isotropic 
ferromagnets. 16-18 ' 23 ' 56 The spin wave stiffness has de- 
creased considerably over this temperature range upon 
approaching Tc- Similar measurements were made to 
determine D in the 18 sample. Figure 10(b) shows 
a comparison of D versus T for the two samples. The 
two samples show identical results for D over the entire 
temperature range where spin waves are observed. Note 
also that there are in fact two complete sets of data for 



the 16 O, one taken with the two samples heat treated in 
parallel, and one set taken after exchanging the 18 O for 
16 O. Both sets of data are identical. The only difference 
found is that for the lower-Tc 18 O sample the measure- 
ments cannot be continued to as high a temperature as 
with the 16 sample. This result obviously cannot be un- 
derstood in the conventional context of fcrromagnetism 
since Tc has changed but the ferromagnetic interaction 
as measured by D has not. Hence it is clear that Tc is 
not determined by magnetic interactions alone. 

The D(T) data can be well described by a power law 
in this temperature range, with an extrapolated value of 
Tc = 271 K which is identical for both isotopes. The ex- 
ponent used in this fit is [3 — v — 0.4, which is fairly close 
to the prediction of the Hciscnbcrg model (3— v = 0.324. 35 
The particular value of the exponent is not the central 
issue here; the point is that the exchange interaction is 
not affected by the isotope substitution, and the extrap- 
olated Tc is identical for both samples. For both sam- 
ples the actual magnetic transition occurs before D has 
renormalized to zero, but note that the renormalization 
has progressed much farther than is the case for samples 
with lower x, lower Tc samples. Indeed, Fernandez-Baca 
et al. 17 have established that D changes little over a wide 
range of tolerance factor, while large changes in Tc oc- 
cur. From this point of view the (O isotope-independent) 
magnetic interactions determine this extrapolated Tc, 
but polaron formation truncates the transition, thereby 
determining the observed Tc- This truncation occurs 
in a first-order manner and leaves D at a finite value, 
and is the determining factor for Tc over a range of x. 
Other measurements of D(T) for smaller x ~ 0.3 com- 
pounds are very similar to the x = 3/8 samples, just 
with a different truncation value at the different Tc's. 
The downturn in Tc with decreasing x in the phase dia- 
gram is then caused by polaron formation rather than a 
weakening of the double exchange mechanism. This be- 
havior continues until the system becomes insulating, and 
D(T = OK) suddenly drops. For our insulating x = 0.2 
sample D(T = 0) is 3 times smaller than in the opti- 
mally doped sample, in good agreement with the results 
of Refs. 37 and 57. 

This rapid truncation of ferromagnetism is clearly re- 
vealed in Fig. 11, which shows the magnetic fluctuation 
spectrum as a function of temperature at a wave vector 
q = 0.11 A -1 for the two samples. The elastic incoher- 
ent nuclear scattering and a flat temperature indepen- 
dent background have been removed from these spectra, 
and the solid curves are fits to the spin wave disper- 
sion [Eq. (1)] plus a quasielastic (spin diffusion) spectral 
weight function convoluted with the instrumental reso- 
lution. Spectra at a series of g's were obtained at each 
temperature, and the data collection time between runs 
at different temperatures is sufficiently long that there 
are no irreversibility effects due to warming and cooling. 
For the 16 O sample [Fig. 11(a)] we see a three component 
spectrum at 265 K, just below Tc- The spin waves are 
observed in energy gain (E<0) and energy loss (E>0) 
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FIG. 11. Constant-g scans at a wave vector of 0.11 A -1 
showing the nature of the magnetic fluctuation spectrum in 
the vicinity of Tc for the ie O and ls O samples, (a) For the 
ls O sample there is a three peak structure at 265 K, while 
just above Tc at 268 there is a single quasielastic peak, (b) 
For the 18 O sample at 258 K most of the weight of the mag- 
netic scattering is in the spin wave excitations. At 261 K 
the weight in the spin wave scattering has decreased, and the 
spectral weight has shifted to the central component. At 
262 K the scattering is completely diffusive in nature. There 
is no evidence of the usual collapse of the spin wave ener- 
gies found in a second order transition. Note that the count 
rates for the two samples are different even though the sam- 
ples were approximately the same size, and this difference is 
due to the data being collected under different experimental 
conditions; 16 data were taken on BT-2 with collimation 
of 20'-10'-10'-20', while the ls O data were collected on BT-9 
with collimation of 10'-13'-10'-16'. 

at ~0.8 meV. At 268 K, just above Tc, the spin waves 
have disappeared and the scattering is purely quasielastic 
(peaked at E = 0). For the 18 O sample [Fig. 11(b)], on 
the other hand, we observe well defined spin waves below 
Tc at 258 K along with a quasielastic component fitted to 
the data. At 261 K the central component has increased 
to form a clear peak (at this q), while the spin waves 
have lost over 50% of their intensity and have moved to 
lower energies, consistent with all previous observations 
of the lower x compounds which show growth of the cen- 
tral peak at the expense of the spin waves. But further 
increasing temperature by only 1 K wipes out the spin 
waves entirely, leaving only a central peak that corre- 
sponds to ferromagnetic correlations in the paramagnetic 
phase. This very rapid temperature evolution is evidence 
for the truncation of the magnetically ordered phase in 
a first-order manner. The nature of the magnetic fluc- 



tuation spectrum is the same for both samples, except 
that the transition for the 16 O sample is extended higher 
in temperature by 5 K, which results in a D(Tc) value 
that is about 30% lower than for the 18 sample. Note 
in particular that the data for the 18 O sample is pure 
spin diffusion for 262 K, while for the 16 sample we 
have propagating excitations at 265 K. We remark that 
we have taken data for temperatures between the ones 
shown in Fig. 11 to examine the possibility that there is 
a sudden decrease in the spin wave energies that might 
be attributable to a sudden decrease in the effective ex- 
change energy, but all the data indicate that we have 
two-phase coexistence, 16 a ferromagnetic-metallic phase 
and a paramagnetic-polaron glass insulator. For lower 
x samples the truncation occurs at lower temperature, 
with a corresponding higher value of D(Tc) as indicated 
in Fig. 10(b). Thus the transition at lower x is more 
strongly first order. We have also found that the central 
peak usually was observable over a wider range of tem- 
perature at smaller x, 15 K or more, compared to the 
~5 K range observed here. This may also be related to 
how close the system is to the extrapolated second-order 
phase transition temperature Tc. 43 ' 56 

C. Phonon Density-of-States 

With the strong evidence that polaron forma- 
tion plays a crucial role in the ferromagnetic- 
metallic/paramagnetic- insulator phase transition, the 
question becomes what controls polaron formation. 
Given the general relationship between isotope effects 
and lattice dynamics, and the electron-phonon coupling 
through the Jahn- Teller effect, 58 we decided to measure 
the phonon density-of-states in our samples using time- 
of-flight and filter analyzer spectrometer techniques. 

The time-of-flight Fermi chopper spectrometer uses 
cold monochromatic neutrons incident on the sample, 
and then detects the inelastically scattered neutrons in 
a bank of detectors. Near room temperature one can re- 
liably measure energy transfers up to ~ 30 meV. The 
density-of-states is obtained by removing the thermal 
population factor, weighting the scattering by q 2 (as- 
sumes phonon scattering), and then averaging over a 
range of q. Note that for low energies and modest values 
of q a considerable part of the spectrum is magnetic in 
origin for these samples. We have made no attempt to 
separate out the magnetic contribution from the lattice 
(nuclear) contribution. 

The filter analyzer spectrometer (FANS) uses cooled 
polycrystalline beryllium and polycrystalline graphite as 
an analyzer, so that only neutrons with energies below 1.7 
meV reach the detectors. The detectors span a large solid 
angle (1.3 steradians), and the incident energy is scanned 
to determine the inelastic scattering. In contrast to the 
Fermi chopper spectrometer this instrument is designed 
to operate in energy loss mode and no correction needs 
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FIG. 12. Phonon density-of-states for (a) ie O and (b) ls O 
samples taken above and below Tc- Units are those appro- 
priate for the FANS spectrometer with a Cu(220) monochro- 
mator. 

to be is made for q 2 since a 1/q 2 term enters into the 
resolution volume. A run with an empty sample can 
was used to remove background. We used a pyrolytic 
graphite monochromator over the range of 15 meV to 40 
meV and a Cu(220) monochromator for energies between 
38 meV and 120 meV. Data sets for each sample were 
merged together after background subtraction. A similar 
merging was performed to include the low-energy time- 
of-flight data. After merging, the data were corrected for 
multiphonon effects, which tended to be significant only 
above 75 meV. Additional experimental details are given 
in Ref. 59. 

Results for the phonon density-of-states taken above 
and below T c appear in Fig. 12. Identifiable peaks ap- 
pear in the data at 17, 25, 42, 51, and 63 meV. There 
are no phonons above 85 meV after multiphonon cor- 
rections. These observations are in agreement with re- 
cent single-crystal measurements of the phonon spectra of 
Lao.7Cao.3Mn03, which show no phonons above 72 meV 
and substantial flat regions in the dispersion branches 
at 23 meV and 50 meV and an optical band near 40 
meV. 60 Corresponding peaks in Raman scattering for an 
x = 0.33 sample are found at 17, 25, and 51 meV. 61 
Other Raman peaks are also observed at 8, 10, 30, and 
54 meV with a broad maximum at 85 meV. Single crys- 
tal phonon dispersion measurements on the related CMR 
compound Lao.7Sro.3Mn03 show phonon branches up to 
72 meV with no gaps between branches when averaging 
over the entire zone. 62 Substantial magnon-phonon in- 
teractions have also been reported, 63 but it is likely that 
these effects can only be observed in single crystals. For 
the present density-of-states data at low energies there 
are substantial differences above and below Tc, but this 
originates from the change in the magnetic scattering 
from spin diffusion above Tc to propagating spin waves 
below. At higher energies the biggest difference is in the 
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FIG. 13. Comparison of the phonon density-of-states for 
16 O and 18 O samples above Tc- There are some slight differ- 
ences which are more easily seen at low temperature, shown 
in the inset. 

16 data between 75 and 80 meV, where the signal 
at 305 K is slightly higher than at 225 K. A similar ef- 
fect can be seen in the 18 O data except that the scatter 
in the data is larger. In other perovskites this energy 
region corresponds to the oxygen breathing modes but 
there was no specific peak seen in our data that could 
be identified with a mode of vibration in this energy re- 
gion. There may be a slight decrease in intensity of the 
peak near 51 meV in the 16 O data. Both these trends 
are borne out by measurements at 10 K in the 40 to 100 
meV range. The published single crystal results 60 for 
the x = 0.3 sample are far more conclusive and show sig- 
nificant damping effects on the longitudinal Jahn- Teller 
(100) breathing mode which has an energy of 72 meV at 
the zone center and 51 meV at the zone boundary. Near 
the zone boundary the intensity of the phonons strongly 
decreases with temperature until the mode disappears 
above Tc- This is consistent with a decrease in the 51 
meV peak seen in our data. Similar changes are also seen 
in the Raman scattering. 61 Overall, however, the changes 
in the phonon density-of-states are modest. 

A comparison between the different samples above Tc 
also reveals that they are nearly identical, as shown in 
Fig. 13 where we compare the phonon density-of-states 
between samples rather than between temperatures. The 
different temperatures were chosen to reflect the differ- 
ent Tc values. The region that looks most affected by 
oxygen isotope substitution at these temperatures is be- 
tween 20 and 40 meV where the scattering intensity near 
30 meV is lower for the 18 sample than the 16 sam- 
ple. Raman measurements by Franck et al. 51 near K 
showed that the oxygen modes at 28 meV and 53 meV 
are shifted downward by 1 meV and 2.6 meV for the 18 
sample. At low temperatures (inset) there does seem to 
be a difference in the spectra at 68 meV where the 16 O 
peak may be slightly higher and more clearly defined. 
Overall, however, only small changes are observed in the 
density-of-states, either with change in isotope, or with 



change in temperature above versus below Tc, which is 
in agreement with recent tunneling data. 64 

V. DISCUSSION 



conventional second-order transition, that is uncoupled 
from the resistivity transition. 65,66 It will be interesting 
to determine if there are dynamic polarons in the Sr and 
Ba substituted systems. 



The increase in the ferromagnetic transition when 18 O 
is replaced by 16 O, and the concomitant extension of the 
spin wave regime without any change in the overall mag- 
netic energetics, demonstrates in an unambiguous way 
that the ferromagnetic state is prematurely truncated, 
and that the transition from the ferromagnetic-metallic 
state to the paramagnetic-insulator state is first order 
in nature. For the optimally doped x = 3/8 composi- 
tion the extrapolated Curie point, where the transition 
would be second order, is just above the observed fer- 
romagnetic transition, indicating that the transition is 
weakly first order at this composition. It is this proxim- 
ity to the second-order fixed point that enables the large 
(~30%) decrease in the value of the spin stiffness upon 
18 O replaced by 16 O, making it unambiguous that the 
ferromagnetic state is truncated. At somewhat higher 
Ca concentrations the transition in fact becomes second 
order. 43 

The origin for the first-order truncation is the sudden 
formation of polarons in the system, which trap the car- 
riers and render the system insulating. It is the polaron 
formation temperature that determines the observed Tc 
in the first-order regime, not the strength of the mag- 
netic interactions. Replacing 16 O by 18 O lowers the en- 
ergy of the Jahn- Teller phonon, making it easier to form 
polarons and consequently lowers the transition temper- 
ature. The application of stress also lowers the barrier 
for polaron formation. Lowering the Ca concentration 
below the x = 3/8 composition decreases the observed 
Curie temperature, while the exchange energy is (ini- 
tially) unaffected. 17 This downturn in Tc (observed) in 
the phase diagram then originates from lowering the en- 
ergy for polaron formation, as it becomes easier for the 
system to form correlated Jahn- Teller distortions; recall 
that at x = pure LaMn03 is strongly Jahn- Teller dis- 
torted, forming a structure with long range orbital (and 
antifcrromagnetic) order. The temperature range where 
the polaron glass is observed 26,27 would then be expected 
to become larger with decreasing x. This trend contin- 
ues until the system becomes insulating around x « 0.2, 
where the double-exchange interaction is lost and the 
magnetic energetics drops precipitously. 37,57 

Ultimately it is the structure that determines the en- 
ergetics for polaron formation, which can be cast in the 
form of the tolerance factor t. 40,47 The Ca system has 
a relatively small t and forms polarons relatively eas- 
ily, and the polaron formation occurs in the ferromag- 
netic state and causes a first-order transition. Larger 
ions such as Sr or Ba have a substantially larger t, and 
do not appear to form a significant number of polarons 
below Tq- The ferromagnetic transition then becomes a 
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